Adiabatic loading of a Bose-Einstein condensate in a 3D optical lattice 
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We experimentally investigate the adiabatic loading of a Bose-Einstein condensate into an optical 
lattice potential. The generation of excitations during the ramp is detected by a corresponding 
decrease in the visibility of the interference pattern observed after free expansion of the cloud. 
We focus on the superfluid regime, where we show that the limiting time scale is related to the 
redistribution of atoms across the lattice by single-particle tunneling. 

PACS numbers: 03.75.Lm,03.75.Hh,03.75.Gg 
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The observation of the superfluid to Mott insulator 
(Ml) transition Q undergone by an ultracold Bose gas 
in an optical lattice has triggered a lot of experimen- 
tal and theoretical activities (see H H 0])- This sys- 
tem allows to experimentally produce strongly-correlated 
quantum systems in a well controlled environment, with 
applications ranging from the realization of novel quan- 
tum phases in multi-component systems (see, e.g., |4( and 
references therein) to the implementation of collisional 
quantum gates HH0. 

Most of these proposals rely on producing a system 
very close to its ground state. However, in experiments 
so far, the successful production of an ultracold gas in the 
optical lattice relies on adiabatic transfer. The conden- 
sate is first produced and evaporatively cooled in order to 
minimize the thermal fraction. The almost pure conden- 
sate is subsequently transferred into the lattice potential 
by ramping up the laser intensities as slow as possible 
in order to approach the adiabatic limit. For too fast 
a ramp, excitations are generated in the system, which 
eventually results in heating after the cloud has equili- 
brated at the final lattice depth. 

It is important in this respect to note the existence of 
two energy scales in this problem, related to the single- 
particle band structure on the one hand or to the many- 
body physics within the lowest Bloch band on the other 
hand. Adiabaticity with respect to the band structure 
is associated with the absence of interband transitions. 
The characteristic time scales are on the order of the 
inverse recoil frequency, typically hundreds of microsec- 
onds. Experimentally, adiabaticity with respect to the 
band structure is easily ensured, and can be checked by 
detecting atoms in the higher Bloch bands H . Adia- 
baticity with respect to the many-body dynamics of the 
system involves considerably longer time scales (on the 
order of tens or even hundreds of milliseconds). Theoret- 
ical studies of the loading dynamics have been reported 
in [H El El El El EH It is clear that non-adiabatic 
effects are more pronounced in the superfluid phase. The 
insulator phase is indeed expected to be quite insensitive 
to such effects, due to the presence of an energy gap. In 
this paper, we focus on the superfluid phase. Our goal is 
to to clarify experimentally how slowly the loading has to 



proceed to minimize unwanted excitations and heating. 

To this aim, we make use of long range phase prop- 
erties characteristic to a BEC. A key observable for ul- 
tracold Bose gases in optical lattices is the interference 
pattern observed after releasing the gas from the lat- 
tice and letting it expand for a certain time of flight 
[1 H El El HE El The contrast of this interference 
pattern is close to unity when most atoms belong to the 
condensate, but diminishes with the condensed fraction 
as the cloud temperature increases. Hence, the visibil- 
ity of this interference pattern can be used to investigate 
the dynamical loading of a condensate into the optical 
lattice. We focus here on a specific lattice depth in the 
superfluid regime. We find a characteristic time scale of 
~ 100 ms above which the visibility of the interference 
pattern appears to be stationary. We show how it re- 
lates to the redistribution of atoms in the lattice through 
single-particle tunneling. 

In our experiment, a 87 Rb Bose-Einstein condensate 
is loaded into an optical lattice created by three or- 
thogonal pairs of counter-propagating laser beams (see 
1] for more details). The superposition of the lattice 
beams, derived from a common source at a wavelength 
Xl — 850 nm, results in a simple cubic periodic poten- 
tial with a lattice spacing d = Al/2 = 425 nm. The 
lattice depth Vb is controlled by the laser intensities, and 
is measured here in units of the single-photon recoil en- 
ergy, E r = h 2 /2m\ 2 L w/ix3.2 kHz. The optical lattice 
is ramped up in a time r ramp , using a smooth waveform 
that minimizes sudden changes at both ends of the ramp. 
The ramp form is calculated numerically. The program 
imposes that the lattice depth is initially zero and reaches 
its final value after a time r ramp . In between (for times 
< t < T r amp), it tries to match the functional form 



V (t) = 



1 + exp 



\ Tramp J 



(1) 



with a — 20, by a piecewise linear approximation. The 
function Q has been proposed in |l2( and is shown in 
Fig. [TJi. After this ramp, the cloud is held in the lattice 
potential for a variable hold time ihoid (see Fig. ^,), dur- 
ing which the system can re-thermalize. After switching 
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FIG. 1: (a)Sketch of the time profile used to ramp up the 
lattice depth to its maximum value Vmax. (b) Evolution of 
the visibility of the interference pattern as the lattice depth 
is increased. The set of data shown corresponds to ~ 3 x 10 5 
atoms (grey circles). 



off the optical and magnetic potentials simultaneously 
and allowing for typically t = 10 — 22 ms of free ex- 
pansion, standard absorption imaging of the atom cloud 
yields a two-dimensional map of the density distribution 
n (integrated along the probe line of sight). 

To extract quantitative information from time-of-fiight 
pictures, we use the usual definition of the visibility of 
interference fringes, 



V = 



(2) 



In this work, we follow the method introduced in 0] 
and measure the maximum density n max at the first lat- 
eral peaks of the interference pattern, (i.e. at the cen- 
ter of the second Brillouin zone) . The minimum density 
n m i n is measured along a diagonal with the same dis- 
tance from the central peak (see inset in Fig.^ 3 ). In this 
way, the Wannier envelope is the same for each term and 
cancels out in the division (see ^3)- For a given two- 
dimensional absorption image, four such pairs exist and 
the corresponding values are averaged to yield the visi- 
bility. In addition, it is worth pointing out that this is 
an essentially model-independent characterization of the 
many-body system. 

For comparison, we reproduce here measurements of 
the visibility as a function of lattice depth, similar to 
those presented in 01 • The data corresponds to a given 
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FIG. 2: Time evolution of the visibility for a fixed final depth 
Vmax = 10 E r and various ramp times: (a): T ramp = 160 ms, 

(b): Tramp = 80 HIS, (c): Tramp = 40 HIS, and (d): Tramp = 

20 ms. 



total atom number N w 3 x 10 5 atoms (grey circles in 
Fig. QJd). For lattice depths larger than 12.5 E r , the 
system is in the insulating phase, as demonstrated in 
Yet, the visibility remains finite well above this point. For 
example, at a lattice depth of 15 E r , the contrast is still 
around 30%, reducing to a few percent level only for a 
rather high lattice depth of 30 E r . In [l^, we have shown 
that such a slow loss in visibility is expected in fact even 
in the ground state of the system, being a consequence of 
the admixture of particle/hole pairs to the ground state 
for finite lattice depths. For these experiments, a ramp 
time of 160 ms followed by a hold time of 40 ms were 
used. We will show later on that this is slow enough to 
ensure adiabaticity. 

We now focus on the superfluid regime, at a lattice 
depth Vq = 10 E r . To investigate how the system is af- 
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fected by the ramping procedure, we fix the ramp time 
Tramp and record how the interference pattern evolves as 
a function of hold time Thoid • Examples of such measure- 
ments are shown in Fig. [SJi-d. 

For the slowest ramp shown, r ram p = 160 ms, the visi- 
bility decreases with a time constant ~ 500 ms. In fact, 
a very similar behavior is observed as soon as the ramp 
time exceeds T ramp > 100 ms. Since this behavior is es- 
sentially independent of the ramp speed, it points to the 
presence of heating mechanisms which significantly de- 
grade the visibility on a time scale of several hundred ms. 
The source of heating could be technical, due for instance 
to intensity or pointing fluctuations of the lattice beams, 
or intrinsic processes. An unavoidable process is for in- 
stance atomic spontaneous emission following excitations 
by one of the lattice beams. Here we attempt a crude es- 
timate of the effect of spontaneous emission on visibility 
as follows. The heating rate Thcat ~ r sp x (/i 2 /2toAq) is 
given by the rate T sp at which such events happen due 
to all lattice beams, times the recoil energy 
where Ao ~ 780 nm is the resonant wavelength of the D2 
transition. To estimate the time scale Ty over which the 
visibility vanishes, we compute the time over which the 
energy gain per particle, FheatTy, is on the order of the 
critical temperature T c times Boltzmann's constant ks- 
In a lattice potential with approximately one atom per 
site, we have ksTc ~ zJ, where J is the tunneling ma- 
trix element and z = 6 the number of nearest neighbors 
in a three dimensional cubic lattice. This yields a final 
estimate 

(3) 

For a lattice depth of 10 E r , we calculate a tunnel- 
ing energy zJ/E r ~ 0.12 and a total scattering rate 
r sp ~ 0.2 s , giving Ty ~ 0.6 s, in good agreement with 
the observed time scale. We conclude that on the time 
scales considered here, photon scattering contributes sig- 
nificantly to the observed heating rate. Although heating 
due to technical noise is not excluded, the bound above 
implies that it is not much more severe than photon scat- 
tering. We note that in principle, the decay in visibility 
could be used to measure the heating rate, provided the 
dependance of the visibility on temperature is known. 

When the ramp times is decreased below T ramp = 
100 ms, the visibility decreases in a qualitatively differ- 
ent way. As shown in Fig. |2>d, the decay occurs with a 
much shorter time constant. Note that at long hold times 
Thold = 800 ms, the visibility has dropped to V ~ 0.2, 
a value almost independent of the ramp time. This is 
consistent with our earlier claim that heating dominates 
at long times. We attribute the short-times decay to the 
generation of excitations by a finite ramp speed. To com- 
pare the different ramp times, we plot in Fig. the mea- 
sured visibility as a function of ramp time for a fixed hold 
time of 300 ms, long enough for the excitations generated 
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FIG. 3: Measured visibility versus ramp time, for a fixed 300 
ms hold time. 
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FIG. 4: Measured visibility decay time versus ramp time. 
The time constant is the half-width at half-maximum of a 
Lorentzian fit to the data. Error bars indicate the statistical 
error. The solid line is a fit to an inverted Gaussian, returning 
a time scale of ~ 60 ms for the visibility decay time to become 
independent of the ramp. 

during the ramp to relax, but short enough that heating 
effects do not blur all differences between the different 
ramps. A time scale of ~ 100 ms clearly emerges, above 
which the ramp time can be increased without notice- 
able effect. To check whether this could be an artifact 
of our plotting method, we have fitted the data to an 
empirically chosen Lorentzian form. The half-width at 
half-maximum of the Lorentzian was taken as the visibil- 
ity decay time. As seen from Fig. 0] the same behavior 
is found, with a characteristic 1/e ramp time w 60 ms. 

This time scale can be understood using elementary 
arguments. Adiabatic evolution in a quantum system 
with time-dependent hamiltonian requires the condition 

\H\ « h\u> fi \ 2 , (4) 

where H is the derivative of the hamiltonian and where 
LUfi is the Bohr frequency of transition between the (in- 
stantaneous) eigenstates |z) and |/). Eq. J3J should be 
fulfilled at all times. For an ultracold gas in an optical 
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lattice, three energy scales appear: (i) the tunneling ma- 
trix element J defining the rate of hopping from site to 
site, (ii) the on-site interaction energy U and, (iii) the en- 
ergy associated to the "external" confinement potential 
usually present on top of the lattice. This potential V ex t 
results from the combination of the magnetic trapping 
potential in which the condensate is initially produced, 
and of an additional confinement due to the Gaussian 
shape of the lasers producing the optical lattice. In prac- 
tice, V ex t is nearly harmonic, with a trapping frequency 
SI = u> cxt ~ y/^m + 8Vo/mw 2 , where U! m = 2ir x 16 Hz is 
the oscillation frequency in the magnetic trap and where 
w w 136 /im is the laser beam size. With a proper choice 
of the laser beam sizes, the increase in U and f2 are such 
that the equilibrium size of the condensate as calculated 
in the Thomas-Fermi approximation varies little during 
the ramp |20j . In addition, when the lattice depth is 
increased, J drops exponentially fast whereas U and f2 
increase slowly. Hence, the question of adiabaticity es- 
sentially reduces to whether the atoms can redistribute 
through tunneling in order to adapt the size of the system 
to the instantaneous Thomas-Fermi shape. This suggests 
to take \H\ ~ \ J\ and uj/i ~ J/h in Eq. @, giving the 
following definition for the adiabaticity parameter A, 



A 



Max 

0<t<T r ,,„ 



h\j\ 
I 2 



(5) 



For a given ramp, we require A <C 1 to ensure an adia- 
batic loading of the cloud into the lattice. In the inset 
of Fig. we show the quantity h\J\/J 2 calculated for 
50 ms. We have calculated this curve for a fi- 



' ramp 



nal lattice depth of IOEr using the lattice ramp given in 
Eq. m, and the analytical estimate 0, 



J 



^{Vo/E r f^ exp{-2^V /E r ). 



(6) 



The sharp decrease for short times in the inset of Fig. 
follows from the inadequacy of the tight-binding approx- 
imation under which Eq. © is derived. We ignore this 
feature, and calculate A from the peak value occurring 
near t « r ramp /2, where the rate of change of the lattice 
depth is highest. We have repeated the calculation for 
several ramp times (see Fig- - We find that the critical 
A = I corresponds to a ramp time r ramp « 80 ms, close 
to the experimental findings. This suggests that for our 
experimental parameters, the loading is indeed limited 
by single-particle tunneling. 

In conclusion, we have studied how the visibility of the 
interference pattern was affected by the speed at which 
the lattice was ramped up. A time scale of ~ 100 ms 
was found for adiabatic loading in the optical lattice . 
In this paper, we focused on the dynamics properties 
in the superfluid regime. Even more interesting is the 
dynamical evolution of the system as the superfhiid-to- 
Mott-insulator transition is crossed. An important and 
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FIG. 5: Adiabaticity parameter (see text) plotted versus ramp 
time. In the calculation, we have used the actual ramp form 
as generated in the experiment to describe the lattice depth 
increase. In the inset, we show how the adiabaticity parame- 
ter changes in time as the lattice depth is ramped up. A ramp 
time of T ramp = 50 ms has been used for this plot. 



still open question is in particular how reversible this 
transition is. Experiments 0, 0, H^ ? ] found that 
one could ramp up the lattice intensity and reach the 
regime where phase coherence is lost, then ramp down 
the lattice and regain it back. To what extent the ini- 
tial phase coherent state can be recovered, and what are 
the limiting mechanisms has however not been studied. 
Also, the recent availability of virtually exact methods 
for one-dimensional systems |l2j could allow for a pre- 
cise comparison with the experimental data, in a strongly 
non-equilibrium regime where theory is still in progres- 
sion. Finally, new dynamical effects are predicted, such 
as oscillation of the order parameter and vortices-driven 
relaxation analog to the Kibble-Zurek mechanism |2l| . 
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